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Abstract 

We investigate the magnetotransport properties of quasi-free standing epitaxial graphene bilayer 
on SiC, grown by atmospheric pressure graphitization in Ar, followed by H2 intercalation. At the 
charge neutrality point the longitudinal resistance shows an insulating behavior, which follows a 
temperature dependence consistent with variable range hopping transport in a gapped state. In 
a perpendicular magnetic field, we observe quantum Hall states (QHSs) both at filling factors (y) 
multiple of four [y — 4, 8, 12), as well as broken valley symmetry QHSs at v — and v — 6. These 
results unambiguously show that the quasi-free standing graphene bilayer grown on the Si-face of 
SiC exhibits Bernal stacking. 
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Graphene bilayers in Bernal stacking; 111 exhibit a transverse electric field tunable band- 

nn 

gap \Zu lol l , as evidenced by angle-resolved photoemission |4| and transport measurements 
P,l6|, a property that renders this material attractive for device applications. Bernal stacking 
is the lowest energy structure, and is found in the natural graphite crystal. Recent studies 
of graphene bilayer grown on SiC [3] and Cu [8|, |9|] have suggested the presence of Bernal 
stacking, based primarily on electron microscopy and Raman spectroscopy. However, there 
has been no firm evidence of Bernal stacking based on electron transport, which in turn 
would allow an assessment of its potential role for electronic devices [10j. It is therefore 
of high interest to explore the electron transport in grown graphene bilayers, in order to 
determine the stacking of the layers and the key transport properties. 

Here we investigate the transport properties of epitaxial graphene bilayer grown by at- 
mospheric pressure graphitization of SiC followed by H 2 intercalation, which renders the 
graphene quasi-free standing. Using top gated Hall bars with AI2O3 dielectric we probe 
the magneto-transport up to magnetic fields of 30 T and temperatures down to 0.3 K. The 
devices show a high field-effect mobility of 2,600 - 4,400 cm 2 /Vs, which changes little from 
room temperature down to 0.3 K, as well as a strong insulating behavior near the charge 
neutrality point. The magneto-transport data reveal quantum Hall states (QHSs) at filling 
factors v — 4, 8, 12, consistent with the four- fold, spin and valley degenerate Landau levels 
(LL) in Bernal stacked (A-B) graphene bilayer. More interestingly, the data also reveal 
developing broken valley symmetry QHSs at filling factors v — and v = 6, which testify to 
the high sample quality. Further supported by Raman spectroscopy and low energy electron 
microscopy (LEEM) data, these results unambiguously show that the quasi-free-standing 
epitaxial graphene bilayer grown on the Si face of SiC substrates exhibits Bernal stacking. 

The graphene bilayer films studied in this paper are produced via a two step process 
beginning with a starting substrate of 6H-SiC(0001) (Si-face, 2.1 xlO 11 Q-cm, II- VI Incor- 
porated). Prior to graphitization, the substrate is hydrogen etched (45% H 2 - Ar mixture) 
at 1350 °C to produce well-ordered atomic terraces of SiC. Subsequently, the SiC sample is 
heated to 1000 °C in a 10 % H 2 - Ar mixture, and then further heated to 1550 °C in an 
Ar atmosphere {ll, 12]. This graphitization process results in the growth of an electrically 



active graphene layer on top of the buffer layer, covalently bound to the substrate. Finally, 
hydrogen intercalation was carried out using 45% H 2 - Ar mixture at 800 °C [7j, in order 
to decouple the buffer layer from the substrate. As we show here, the two graphene layers, 
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a formally buffer layer (decoupled via hydrogen intercalation) and a monolayer graphene 
(formed via graphitization process), are in Bernal stacking. 

Both the number of layers and the quality of the resulting material are probed by Raman 
spectroscopy acquired over a 25x25 fim region, using a 532 nm excitation wavelength, 5 
mW power, and 500 nm spot size. Figure 1(a) shows a typical Raman spectrum of the 
graphene bilayer sample in which the contribution from SiC substrate has been subtracted. 
Well defined spectral features characteristic of graphene's G (^1600 cm -1 ) and 2D-band 
(^2700 cm -1 ) are observed. The 2D-band was well fit utilizing four Lorentzian peaks each 
with a width of ^35 cm -1 confirming that the presence of graphene bilayer [Fig. 1(a) inset]. 
The Raman spectra indicate that the sample is of high quality as the intensity ratio of the 
defect induced D-band to G-band, while present is less than 0.1 for the great majority of 
the sample area [Fig. 1(b)]. The presence of graphene bilayer is further confirmed by the 
energy dependence spectrum of the specular electron reflection [l^t , and by the well ordered 
6- fold low energy electron diffraction (LEED) pattern acquired with LEEM [Fig. 1(c)]. No 
turbostatic disorder was observed. 

To probe the transport properties of these graphene bilayers we fabricate top-gated Hall 
bars. The Hall bar location on the substrate is first chosen using optical and atomic force 
microscopy (AFM) in order to identify an appropriately wide terrace [Fig. 1(d)]. The 
dashed contour of Fig. 1(d) is an example of an 8 /xm-wide terrace, onto which a Hall bar 
is subsequently defined. Raman imaging confirms that there is only graphene bilayer within 
these terraces, as panel (b) inset shows that the 2D mode is rather uniform inside terraces. 
E-beam lithography and O2 plasma etching are used to pattern the Hall bar active area; the 
graphene is etched outside the Hall bar to prevent parallel current flow. Metal contacts are 
realized using a second e-beam lithography step, followed by a 40 nm Ni deposition and lift- 
off. To deposit the AI2O3 gate dielectric, we first deposit a 1.5 nm thick Al film. The sample 
is then exposed to ambient, and transferred to an atomic layer deposition (ALD) chamber. 



The ambient exposure causes the Al interfacial layer to fully oxidize 14], and provides 
nucleation centers for the subsequent ALD process. A 15 nm thick AI2O3 top dielectric film 
is then deposited using trimethylaluminum as Al source, and H 2 as an oxidizer [j^j]. A 
third e-beam lithography step, followed by Ni deposition and lift-off are used to pattern the 
top gate [Fig. 2(a) inset]. The corresponding dielectric capacitance for this stack is 245 
nF/cm 2 , with an average dielectric constant of k — 4.6. 
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Four point longitudinal (p xx ) and Hall (p xy ) resistivity measurements are performed at 
temperatures (T) down to 0.3 K, and magnetic fields (B) up to 30 T, using low-current, 
low- frequency lock-in techniques. The carrier density (n), and its dependence on the top 
gate voltage (Vtg) are determined from Hall measurements, as well as the filling factors of 
the quantum Hall states observed in high magnetic fields. The longitudinal (a xx ) and Hall 
(cr X y) conductivities are determined via a tensor inversion from the measured resistivities. 

In Fig. 2(a) we show a xx vs. Vtg^ measured at T = 290 K and 0.6 K, and at B = T, 
revealing an ambipolar characteristic with a charge neutrality point at a positive Vtg value. 
Away from the charge neutrality point, the a xx vs. Vtg data show a linear dependence up to 
the highes t Vt g values, with a corresponding field-effect mobility of 2,600 - 4,400 cm 2 /Vs at 
T=290 K la]. At T = 0.6 K, where the measurement was performed in a wider Vtg range, 
the a xx vs. Vtg linear dependence persists down to -4 V, for a gate voltage overdrive of up 
to 5 V. The linear a xx vs. Vtg dependence contrasts data reported in mono-layer graphene, 
where neutral impurity scattering, which is density independent, limits the conductivity at 
high Vtg values. In contrast, the neutral impurity scattering remains proportional to n in 
graphene bilayers [171 ]. 

Figure 2(b) shows p xx vs. Vtg a t different T values, measured at B = T. Remarkably, 
p xx measured near the charge neutrality point is strongly temperature dependent, with an 
insulating behavior. The insulating phase at the charge neutrality point suggests an energy 
gap, and contrasts data from mono-layer graphene with comparable mobility, where p xx at 
the charge neutrality point is weakly temperature dependent. In contrast, graphene bilayers 
in the presence of an applied transverse electric field open a tunable energy gap between 
the conduction and valence bands, thanks to the layer on-site energy asymmetry [21]. We 
posit the presence of a transverse electric field in our samples, due to unintentional doping 
and the asymmetry of the device structure. The inset of Fig. 2(b) shows the T-dependence 
of the resistivity measured at the charge neutrality point (pnp)- The data clearly follow a 
oc e ( T o/ T ) 1/3 dependence for T lower than 100 K, indicating that variable range hopping rather 
than thermally activated conduction controls the electron transport at low temperatures. 
This has been attributed to disorder-induced localized states in the gap, which reduce the 
effective energy gap and render the T-dependence of the p^p weaker than the 

exponential oc e A / 2/C£?T dependence expected for a band insulator with an energy gap A; ks 
is the Boltzmann constant. The extracted T value corresponding to Fig. 2(b) inset data is 
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0.6 K, similar to previously reported values on exfoliated graphene bilayers Q. A fit of 
the Pnp vs. T data of Fig. 2(b) inset to the exponential oc e A / 2/C£?T dependence for T > 120 
K, yields an energy gap at the charge neutrality point of A = 20 ± 6 meV. 

While Fig. 2 data at the charge neutrality point show an energy gap, suggestive of a 
graphene bilayer, the most important finding of this study is presented in Fig. 3. Figure 3(a) 
shows p xx and p xy vs. Vtg, measured at a high magnetic field of B = 30 T, and a temperature 
T = 0.3 K. The data reveal strong QHSs at integer filling factors v — 4, and v = 8, marked 
by vanishing or local minima in p xxi and corresponding p xy plateaus. The QHSs at integer 
filling factors multiple of four in Fig. 3 stem from the four-fold degeneracy associated with 
spin and valley of each Landau level, and unambiguously identify the material as being a 
Bernal stacked (A-B) graphene bilayer. In contrast, thanks to the 7r-Berry phase of chiral 
Dirac fermions, the strongest QHSs in mono-layer graphene are present at integer fillings 



v = 2, 6, 10, etc.. 



221 ] Figure 3(b) data show a xy vs. Vtg at different values of the 
S-field. The a xy plateaus at v • e 2 /h values, which become stronger with increasing the 
£>-field confirm the presence of QHSs at fillings v — 4, 8, 12 consistent with a Bernal stacked 
graphene bilayer. More interestingly, Fig. 3 data reveal a strong QHS at v — 0, and a 
developing QHS at v = 6. While a possible explanation for the v — 6 QHS is a mixture of 
graphene mono- and bilayer, this is ruled out by the absence of the other QHSs associated 
with mono-layer graphene, in particular the v — 2 QHS. The absence of graphene monolayers 
is further corroborated by LEEM measurements. On the other hand, the presence of v — 
and v = 6 QHSs is fully consistent with a graphene bilayer in the presence of a transverse 
electric field, a finding which is also consistent with Fig. 2 data. We expand below on this 
argument. 

At zero transverse electric field (£"), an eight-fold degenerate Landau level (LL) is located 
at zero energy (charge neutrality point); it consists of the N — 0, 1 LLs with their respective 
valley and spin degeneracies. The N = 0, 1 LLs are layer polarized and split in the presence 
of a transverse E'-field depending on the on-site energy of the layer (valley) degree of freedom 
leading to a QHS at v — 0. Similarly, the presence of a transverse E'-field breaks the 
valley degeneracy of \N\ > 2 LLs [231 ]. but the corresponding energy splitting is smaller 
than that of N = 0, 1 LLs. This explains the absence of v — 6 QHS in exfoliated graphene 
bilayers on Si02, where the LL disorder broadening can obscure the E'-field induced splitting. 
In that regard, the observation of a v = 6 QHS in quasi-free standing graphene bilayers is 
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interesting, and testifies to a reduced disorder and LL broadening in these samples, by 
comparison to dual-gated graphene bilayers on Si02 substrates, and with a similar top-gate 
stack [24j]. We note that back-gated graphene bilayers on Si02 samples can exhibit higher 



mobilities 

In the remainder, we explore further the v = QHS in these samples. Depending on 
the transverse E'-field, the v — QHS in graphene bilayers can be either spin-polarized at 
small E'-fields, or valley- (layer) polarized at large E'-fields (26, 27]. If the E'-field is varied 
at a given perpendicular magnetic field, the v = QHS undergoes a transition from spin-to- 



24|, [28| on dual-gated A-B 



valley-polarized at a critical E'-field (E c ). Experimental studies 
(Bernal) graphene bilayers exfoliated from natural graphite show a linear dependence of E c 
on B ) with a slope of 12-18 (mV/nm)-T _1 . Figure 4(a) shows p xx vs. Vtg, measured at B 
= 30 T, and at different temperatures. The data reveal an insulating phase near the charge 
neutrality point, consistent with a LL splitting at zero energy responsible for the presence 
of the v = QHS of Fig. 3(b). The inset of Fig. 4(a) shows the Arrhenius plot of p xx 
vs. T _1 at v — 4 QHS; the data follow an activated T-dependence, p xx oc e ~^/ 2k BT w ith 
an energy gap A 4 = 16 K for v = 4. Figure 4(b) shows p xx vs. Vtg measured at T = 0.3 
K and at different values of the £>-field. Figure 4(b) data show a weak dependence of p xx 
on B at v = 0. This finding, analyzed by comparison with existing data in exfoliated A-B 
bilayers [24] suggest that the v = is neither spin-polarized or valley-polarized in the range 
of S-fields probed here, but rather in the vicinity of the spin-to- valley-polarized transition 
(see supplementary material). Indeed, if the v = QHS was spin polarized, we would 
expect p xx at v — to increase with B. Conversely, if the v = QHS was valley polarized 
a decrease of p xx with the £>-field is expected. This argument provides us with an estimate 
of the transverse i?-field across the bilayer of 0.35±0.07 V/nm. 

In summary, we investigate the magnetotransport in quasi- free-standing graphene bilayers 
on SiC. We observe QHSs at fillings v = 0,4,6,8,12, consistent with a Bernal stacked 
graphene bilayer in the presence of a transverse field. An insulating state observed at B = 
T, near the charge neutrality point indicates the opening of an energy gap, in agreement with 
the expected response of a Bernal stacked graphene bilayer. These findings corroborated 
with earlier microscopy studies j3] unambiguously identify the Bernal stacking arrangement 
of graphene bilayers on the Si-face of SiC substrates, and render this system particularly 
attractive for electronic and optoelectronic device applications, thanks to its high mobility, 
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tunable energy gap, and high on/off ratio. 
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FIG. 1: (color online) Raman, LEED and AFM characterization of the graphene bilayer sample, 
(a) Representative spectrum of the bilayer region, with the SiC response removed. The inset 
shows a 2D band fitted well using four Lorentzian functions, an indication of the graphene bilayer 
presence, (b) Histogram of the 1(D) /1(G) ratio acquired from Raman mapping (25x25 /mi 2 , 75x75 
data points). The inset shows the spatial distribution of the total 2D peak width indicating that 
the bilayer is present on the terraces (dark regions within inset, total width cm -1 ), whereas 
thicker graphene layers (bright regions of inset, total width ~70 cm -1 ) are located at atomic steps 
originating from the SiC. (c) LEED pattern of a bilayer graphene obtained at the illuminating 
electron energy of 48 eV. (d) AFM topography of the graphene sample. The dashed contour on 
the plateau indicates the region used for device fabrication. 
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FIG. 2: (color online) (a) a xx vs. Vtg measured at T = 290 K and T = 0.6 K. The top axis 
represents the carrier density (n); positive values correspond to n-type carriers (electrons), while 
negative values correspond to p-type carriers (holes). The inset shows an optical micrograph of 
the top-gated Hall bar. (b) p xx vs. Vtg measured at different T values. The resistivity near the 
charge neutrality point shows a clear insulating behavior. The inset shows pisip vs. T -1 / 3 on a 
log-lin scale. The data are consistent with variable range hopping transport, measured in exfoliated 
bilayers when a band gap opens as a result of an applied transverse electric field. 
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FIG. 3: (color online) (a) p xx and p xy vs. Vtg (bottom axis), and n (top axis), measured at B = 30 
T, and T = 0.3 K. (b) a xy vs. Vtg (bottom axis) and n (top axis) measured at T = 0.3 K, and 
different £>-field values. The data show developing QHSs with increasing the 5-field. 



11 



(a) n(10 12 cm" 2 ) 




(b) n(10 12 cm 2 ) 





■10 -8 -6 -4 -2 2 4 6 


50 


i i i i i i i i i i i i i i i i i i 
- T=0.3K i 

. — 15T a 


40 


20 T 1 


30 


i 


20 




10 



- i . i i l l^T I i " 



V T G ( V ) 

FIG. 4: (color online) (a) p xx vs. Vtg (bottom axis) and n (top axis), measured for different T 
values, at B = 30 T. Inset: p xx vs. T -1 at v = 4 on a log-lin scale. The i/ = 4 energy gap is 
A4 = 16 K. (b) p xx vs. Vtg a ^ different S-field values measured at T = 0.3 K. Concomitantly 
with developing QHSs at increasingly higher £>-field values, the p xx measured at v — changes 
little with the £>-field, suggesting the v — is not strongly spin or valley polarized in the range of 
£>-fields explored here. 
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